Singularity in the boundary resistance between superfluid 4 He and a solid surface 
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We report new measurements in four cells of the thermal boundary resistance R between copper 
and 4 He below but near the superfluid-transition temperature T\. For 10~ 7 < t = 1 — T/T\ < 10 -4 
fits of R — R(,t Xb + Bo to the data yielded xt ~ 0.18, whereas a fit to theoretical values based on 
the renormalization-group theory yielded Xb = 0.23. Alternatively, a good fit of the theory to the 
data could be obtained if the amplitude of the prediction was reduced by a factor close to two. The 
results raise the question whether the boundary conditions used in the theory should be modified. 
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An important, as yet unresolved, issue in condensed- 
matter physics is the nature of the boundary conditions 
which are appropriate to describe the interface between 
a liquid and a solid. Experimentally it is difficult to in- 
vestigate this problem because the boundary layer usu- 
ally is of atomic dimensions and contributes negligibly to 
macroscopic properties. This situation is changed near 
a critical point, where the thickness of a boundary layer 
is of the same order as the correlation length £; the cor- 
relation length diverges at the critical temperature. A 
system particularly suitable for the investigation of the 
boundary layer is the superfluid transition of He. 

In bulk superfluid 4 He heat transport occurs through 
the counterflow of the normalfluid (j n ) and superfluid 
(j s ) currents and does not produce a temperature gra- 
dient [Q. However, Landau predicted ffl that a heat 
current Q orthogonal to a solid wall should produce a 
boundary layer in the fluid adjacent to the wall with a 
temperature difference ATj across it. The reason for this 
is the suppression of the order parameter, and thus of 
j s and j n , near the wall. The phenomenon is illustrated 
schematically in Fig. [j]. Within this boundary layer, sup- 
pression of the counterflow implies that some of the heat 
must be carried by diffusive processes. Consequently, a 
thermal gradient is developed within it. However, the 
thermal resistance Rb — ATb/Q is unobservably small 
deep in the superfluid phase where the boundary-layer 
thickness is of atomic dimensions. Only very close to the 
superfluid transition temperature T\ where t = 1 — T/T\ 
becomes small and £ = £,ot~ L ' becomes large does Rb be- 
come measurable in high-precision experiments. Thus 
this phenomenon was found experimentally only rela- 
tively recently Roughly speaking one can assume 
that the temperature gradient in the boundary layer de- 
cays exponentially from —Q/\ at the wall (where there 
is presumed to be no counterflow at all) to zero deep in 
the superfluid (where all the heat is carried by counter- 
flow) with a characteristic length equal to £, as illustrated 



in Fig. Here A ~ \at~ Xx , with an effective exponent 
x\ ~ 0.43, is the diffusive thermal conductivity of the 
fluid j5|. This leads to the estimate 

R b ~ £/A oc r x » (1) 

where Xb — v — x\ ~ 0.24 A renormalization-group- 
theoretical calculation of Rb carried out by Frank and 
Dohm (FD) @,| agrees with this qualitative picture and 
is expected to give the behavior of Rb(t) quantitatively. 
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FIG. 1. Schematic diagram (not to scale) of the tempera- 
ture variation in the cell. The discontinuity at z = {z = 1) 
is due to the classical Kapitza resistance Rk at the top (bot- 
tom) of the cell, and ATb (which develops over a length scale 
£) is due to Rb- The directions of the currents Q,jn,js, and 
of gravity g are indicated. 

The singular boundary resistance is particularly inter- 
esting because it is a direct consequence of the bound- 
ary conditions (BC) which are relevant at the solid-liquid 
interface. For instance, theoretically popular (but physi- 
cally unrealistic) BC in the study of finite-size effects near 
critical points are periodic ones (PBC) 0. For PBC, 
there is no contribution in the fluid to the boundary re- 
sistance at the solid-liquid interface because the order 
parameter is not suppressed. More realistic BC are of 
the Dirichlet type (DBC). In that case it is assumed that 
the order parameter vanishes at the solid wall. This is 
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the case considered by FD and illustrated in Fig. 0. In 
principle, one can imagine a continuous family of BC flio| , 
and an important question which goes well beyond the 
problem of critical phenomena is which BC are relevant 
to real physical systems. 

We find that excellent fits to new data for four separate 
cells (labeled I to IV) can be obtained with the power- 
law Eq. [j} but that the resulting values of Xb are close to 
0.18 and thus significantly smaller than the theoretically- 
expected value xl h ~ 0.24. Alternatively, a reasonable 
fit of the function predicted by the theory || to the data 
can be obtained if this function is multiplied by a factor 
which is significantly less than one and, for three of the 
cells, close to 0.53. Although a detailed theoretical expla- 
nation of this problem may be required, we believe that 
this discrepancy may be indicative of a need to consider 
modifications of the DBC which have been assumed in 
the theory not only of Rb, but also of finite-size effects 
on the specific heat jll] . 

The apparatus was similar to one described previously 
||,[l2|. We used four different sample cells. Cell I was de- 
scribed elsewhere jl3|,[l4]]. It had a diameter of 1.27 cm. 
The 0.12 mm thick stainless steel sidewall was brazed to 
a top copper piece, leaving no gap between it and the 
copper Jl5| . The copper surface at the sample top was 
finely machined and had a measured root-mean-square 
(rms) roughness of 1 /im. The fill capillary entered the 
cell bottom. Two sideplanes (SPs) similar to those first 
used by Duncan et al. ||[| were used to measure the 
helium temperature. They pressed around the entire cir- 
cumference against the outside surface of the sidewall. 
Thermometers on the SPs could sense the helium tem- 
perature, unencumbered by any temperature drop across 
a solid-liquid interface. Cells II and III had a similar 
construction. However, the cell bottom was sealed to the 
bottom copper piece with epoxy (brazing) for cell II (III) , 
and the fill lines entered the cells from the top. The bot- 
tom surfaces were polished copper and the rms surface 
roughness was 0.01 (0.05) /im for cell II (III). Cell IV was 
similar to one used in the microgravity project DYNAMX 
fl6f . It had copper top and bottom cylinders of diameter 
2.28 cm which were epoxied into a stainless steel wall of 
thickness 0.075 mm. The two copper surfaces were pol- 
ished and had a rms roughness of about 0.02 /im. The 
wall was penetrated by two copper SPs which made con- 
tact all around the circumference with the helium. The 
thickness of each SP was 0.075 mm. The fill capillary 
entered the cell top. 

The temperatures of the top (T t ), bottom (Tb), and the 
two sideplanes (T st and T s b) were measured with nano-K 
resolution p7| . In each case the fill capillary came from 
an overflow volume containing the liquid-vapor interface 
on an isothermal shield which was kept at a temperature 
near 2.2 K and steady to ±l/zK pTSj ] . The surface not 
penetrated by the fill capillary was used for the measure- 
ments. The sample was 4 He with less than 1 ppb 3 He. 



The calculation of Frank and Dohm || gives the dif- 
(t) - i?£' l (i ) and thus the values of 



ference SRf 1 



Rl h < 



are determined only to within an additive constant. 
The theoretical result has the form 



2/,, 



5Rl h 



dt' £ 
T2X 



K{f,w,u,i) 



(2) 



Here A and £ are the thermal conductivity and correla- 
tion length above T\, and if is a function of the static 
and dynamic coupling constants 7, u, /, and w of Model 
F (l^] for the dynamics of the superfluid transition. The 
coupling constants depend on t and are known from inde- 
pendent fits of the theor y |20[ to the thermal conductivity 
pl| and heat capacity J2l| , [22| ] near T\. Using their tab- 
ulated values p3[ , we evaluated SRf 1 numerically from 
the equations given in Ref. H and to = 0.05. 
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FIG. 2. The resistance R of (from bottom to top) cell I 
(moved up by 0.1 cm 2 K/W), II, III, and IV (moved down 
by 0.2 cm 2 K/W). For cell I the open circles, solid circles, 
open squares, and solid squares are for data taken at 2, 20, 
30, and 40 /xW/cm 2 respectively. The solid lines are the the- 
oretical result, shifted vertically so as to agree with the data 
near t — 10 -4 . The dotted lines are the theory multiplied by 
factors of 0.76, 0.59, 0.51, and 0.52 for cell I to IV respectively 
and shifted vertically. 

In Fig. I the data for R = R b + R B are shown for 
all four cells on a linear vertical scale. For cell I there 
are data for several power densities Q which show that, 
within the experimental resolution, the results do not de- 
pend on Q Q. There is a large sample-dependent back- 
ground contribution Rb = Rk + Rcu due to the usual 
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Kapitza resistance Rk and the series resistance Rc-a in 
the copper parts of the cell (see Fig. |l|). The theoreti- 
cal prediction || of after vertical shifts which cause 
agreement with the data near t = 10~ 4 , is given by the 
solid lines in the figure. Obviously the agreement is not 
good. 
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FIG. 3. (a): Values of Xb derived from fitting Eq. |3| to the 
data for cells I to IV. The dotted line is the theoretical value, 
and the dashed line is an average value based on the exper- 
iments, (b) : Values of Co derived from fitting Eq. ^| to the 
data for cells I to IV. The theoretical value corresponds to the 
top of the graph at Co = 1. The dashed line is an average 
value for cells II to IV. 

Equation [l] suggests that it should be possible to fit 
the data near T\ with a powerlaw. Thus, as an alterna- 
tive interpretation we considered the data for very small 
t where we can represent J?b by a constant. We fitted 
the function 

R = R Q r Xb + R B (3) 

to the data with t < 10 , least-squares adjusting i?o, £&, 
and Rb- For t = 1 — T/T\ we used the zero-current value 
of T\ at the vertical position of the cell surface in ques- 
tion H and f = Tjje + AT b /2. We obtained the values 
of Xh shown in Fig. |^a. The error bars are 95% confi- 
dence limits derived from the random errors and do not 
include possible systemstic errors. An equivalent fit to 
the theoretical data yielded xf 1 = 0.232 which is shown 
as a dotted horizontal line in the figure. We see that the 
uncertainties of x& are largest for cells II and III, and 
that the data for those cells are consistent with the the- 
oretical value. However, the more precise data for cells 
I and IV are clearly inconsistent with the theory. All 
four experimental results are consistent with Xh — 0.185 
which is shown as a dashed horizontal line in the figure. 
For comparison, we mention that a similar fit to the data 
from cell J of Ref. || yielded x b = 0.14 ± 0.54 ||. One 
sees that the uncertainty in these data is large enough 
for them to be consistent with the theory. 
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FIG. 4. Upper panel: The contribution Ri, = R — B\ with 
Bi determined from a fit of Eq. ^] to the data. Open (solid) 
circles are for cell I (IV) with Si = 1.072 (1.610). The solid 
lines are the fits, (b): Deviations of the data from the fits. 

In order to learn more about the difference between 
theory and experiment, we considered the possibility that 
the theory might give the t— dependence of Rb correctly, 
but that the amplitude is in error by a factor of order 
unity. Thus we fitted the function 

R = C SRf + B l (4) 

to the data over the range t < 10~ 4 . Here B\ represents 
the sum of the contributions from the experimental back- 
ground Rb and from the unknown contribution R\ (to) 
in the theory (see Eq. |^). Both Co and B\ were least- 
squares adjusted. We found that the data could be fitted 
reasonably well, but with Co < 1. We obtained the val- 
ues of Co shown in Fig. |^b (the error bars again are 95% 
confidence limits). The fit to the data is quite good, as 
shown by the dotted lines in Fig. |[ The data from cell 
J of Ref. H yielded C = 0.88 ± 0.26 from a similar fit 
p5| . As can be seen, the error bars based on these older 
data easily overlap the theoretical result Co = 1. 

For all four cells, Co is significantly less than unity. The 
results for cells II to IV are consistent with Co — 0.535, 
which is shown by the dashed line in Fig. ||b. This result 
is smaller than expected from estimates of the theoreti- 
cal uncertainty of the amplitude of Rf 1 . We do not know 
why the result for cell I is somewhat higher, but note that 
cell I was the only one which had a relatively rough and 
unpolished surface. To illustrate the quality of the fit to 
the data, we show Rb = R — B\ and CoSRf 1 for cells I 
and IV (which have the highest precision) in Fig. Ha, and 
the deviations of the data from CoSRl h in Fig. [|b. We 
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see that the fit to the data, although quite good, is not 
perfect particularly for cell IV which yielded the most 
precise data. 

We note that the theoretical result was obtained from 
a first-order perturbation calculation. It would be nat- 
ural to look here for the difference between theory and 
experiment. However, the value of K in Eq. |^ is nearly 
independent of t, leaving the exponent of Rb essentially 
unaltered and nearly equal to that of £/A. It seems un- 
likely that this situation would change at higher order. 
The value of K might change by a factor of order one 
and thus affect Co if a higher-order calculation were per- 
formed; but a change by a factor as small as 0.53 seems 
unlikely. This leaves us to suspect the choice of the DBC 
as the cause of the difference between the theory and the 
experiment. It would be very interesting to investigate 
theoretically whether the BC can be modified so as to 
cause the calculation to match the experimental results 
more closely. 
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